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Abstract 

The behaviour of two-loop two-point diagrams at non-zero thresholds corre- 
sponding to two-particle cuts is analyzed. The masses involved in a cut and the 
external momentum are assumed to be small as compared to some of the other 
masses of the diagram. By employing general formulae of asymptotic expansions 
of Feynman diagrams in momenta and masses, we construct an algorithm to de- 
rive analytic approximations to the diagrams. In such a way, we calculate several 
first coefficients of the expansion. Since no conditions on relative values of the small 
masses and the external momentum are imposed, the threshold irregularities are de- 
scribed analytically. Numerical examples, using diagrams occurring in the Standard 
Model, illustrate the convergence of the expansion below the first large threshold. 
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1. Introduction 



In recent years many precision calculations have been performed in the Standard 
Model (SM), which are required to match the impressive experimental accuracy. A good 
overview of the present status has been recently given It is to be expected that even 
more refined evaluations of loop corrections in the SM and extensions thereof will be 
needed in the future. Since one-loop corrections in the SM are well established the focus 
will be on two-loop studies. As physics examples one should mention the quantity Ar 
occurring in the Fermi constant Gp, small angle Bhabha scattering and the Z ^ bb decay 
rate. 

Such calculations are cumbersome in many respects e.g. through the number of dia- 
grams, the involved tensor structures in diagrams and the complexity of basic building 
blocks i.e. scalar two-loop diagram^. 

It is the purpose of this paper to contribute to the techniques for calculation of the 
scalar two-loop diagrams. When all masses in the propagators are non-vanishing there 
exist arguments that they cannot be expressed anymore in terms of known functions 
like polylogarithms. As an example, we can mention the simplest case of a two-loop scalar 
diagram i.e. a self energy with three massive propagators, the so-called "sunset" (alias 
"sunrise" ) diagram, which can be expressed in terms of Lauricella functions or one- fold 
integrals [^, ||, but no simpler results are available so far. 

In such a situation two main strategies can be chosen. One is a numerical approach, 
another an approximative analytical method. Most progress has been made for self- 
energy diagrams. Obviously they are the simplest ones, but they are physics-wise also 
very relevant. 

As to the numerical strategy one still has a number of options. There exist a two- 
dimensional integral representation (see also in |^) originally introduced for the master 
diagram (Fig. la), but later on extended to all (also divergent) two-loop self-energies [^. 
Another way is the use of dispersion integrals [^, |I0|, whereas other methods rely on 
higher dimensional integration techniques [[lT| , §]. 

The analytical approach uses explicit formulae for the asymptotic expansion of Feyn- 
man diagrams in momenta and masses and is based on general theorems on asymptotic 
expansions (see also |TB[ for review). Some recent physical examples of application of 



these formulae can be found e.g. in [|1^. For two-loop self-energy diagrams with general 



masses, expansions in different regions were systematically examined in refs. |TB|, 0, p!7|f|. 

When a few terms of the series give an adequate description this approach will be 
faster than the purely numerical integration method. Moreover, analytic expressions are 
much more convenient to deal with when the values of some masses are not fixed (e.g. 
the Higgs mass) or have rather wide error bars (e.g. the top quark mass). 

The first results of the analytical approach have been obtained outside the particle 
thresholds: a small momentum expansion below the lowest threshold and a large mo- 

^The tensor decomposition of two-loop self-energy diagrams is described in rcf. [H. 

^The small momentum expansion in the three-point case was studied in ref. [ |l8| where also the 
techniques of conformal mapping and Pade approximations were employed to get numerical results beyond 
the threshold(s). 
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mentum expansion above the highest threshold ||T6[Q- The region between the thresholds 
poses a problem since both small and large momentum expansions do not describe the 
behaviour between the lowest and the highest physical thresholds. In a previous paper 
[[T7| it was shown that in the special case of a zero mass threshold the small momentum 
expansion could be extended to the lowest non-vanishing threshold. In this case the ex- 
pansion coefficients involve the zero-threshold cut which appears as powers of ln(— fc^) 
(where k is the external momentum). 

It is the purpose of the present paper to consider cases when one (or two) two-particle 
threshold(s) is (are) small with respect to the other thresholds, but not anymore zero. 
By using asymptotic expansions in the large mass limit it becomes possible to find a 
series converging above the small two-particle threshold. The expansion coefficients now 
contain the two-particle cut(s) associated with the small threshold(s) and therefore the 
non-regular behaviour around the threshold(s) will be exactly described. Thus the analytic 
approach now is substantially extended into regions, which were hitherto inaccessible due 
to the inapplicability of the used series. From numerical examples it can be seen that 
the convergence holds beyond the small particle threshold(s). We note that there are no 
exact results (except for the zero-threshold limit) available for the cases we are going to 
describe here. 

The actual outline of the paper is as follows. In section 2 different small-threshold 
configurations are considered and subgraphs contributing to the asymptotic expansions 
of the corresponding diagrams are described. In section 3 we present some analytical 
results for the lowest terms of the expansions. Section 4 contains a numerical comparison 
between the direct numerical integration of the Feynman diagram and the results of our 
analytic approximations. Conclusions and future prospects are given in section 4. In 
Appendix A we collect some relevant information on one-loop two-point integrals with 
different masses and present some formulae to handle the numerators. In Appendix B we 
present a result for one of the next-to-leading terms of the expansion. 



2. Constructing the expansion 

In Fig. la,b two different topologies of two-loop self-energy diagrams are shown. The 
dimensionally-regularized scalar Feynman integrals corresponding to these diagrams can 
be written as0 

J ({z/,}; {m,}; k) = J J -^—^^^^^^ (1) 

where n = 4 — 2e is the space-time dimension [^], and {DiY^ = {pj — mf + i^Y^ are the 
powers of the denominators (scalar propagators) corresponding to the internal lines of the 
diagrams in Fig. la,b. The momenta Pi flowing through the internal lines are constructed, 
in an obvious way, from the external momentum k and the loop momenta p and q. 

■^For some special cases (when some of the internal lines are massless), the expansions and exact results 
were presented in ^ . 

^In many cases, we adopt the notation used in the previous papers ^ p7[ . 
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We are mainly interested in the cases when all z/'s are integer. As mentioned in 
refs. fl^, |16], |l3 , in this case it is enough to consider the "master" diagram (Fig. la) only, 
because (i) the diagrams with four and three internal lines can be obtained from Fig. la 
by shrinking (reducing to points) some of the internal lines (this corresponds to putting 
the corresponding i/'s equal to zero) and (ii) the diagram in Fig. lb can be reduced to the 
diagrams with four denominators]^. Thus, in what follows we shall imply that the integral 
(|I|) corresponds to the diagram in Fig. la. 

In general, the diagram in Fig. la has two two-particle thresholds (at k"^ = {mi +1714)'^ 
and k"^ = {1712 + rn^Y) and two three-particle thresholds (at k"^ = (mi -|- 777,3 + ^5)"^ 
and k"^ = (7/7,2 + ^3 + ^4)^)- If all these physical thresholds correspond to non-zero 
(positive) values of /c^, the small momentum expansion is nothing but a Taylor expansion 
in k'^. In [|l^ a general algorithm for calculating the coefficients of this expansion was 
developed. Otherwise, if we have one (or more) physical threshold(s) at k"^ = 0, the 
corresponding "zero-threshold" expansion (as n — 4) involves also the terms with ln(— fc^) 
(and even ln^(— /c^) when two physical thresholds vanish). To calculate analytically the 
coefficients of zero-threshold expansion of the diagram of Fig. la, in general theorems 



on asymptotic expansions of Feynman diagrams |12[ were applied. In such a way, all 
different configurations of massless thresholds have been considered. In particular, it has 
been verified that the expansion converges up to the first non-zero threshold. 

In this paper, we study the behaviour of two-loop self-energy diagrams which have 
small (but non-zero) physical thresholds. Namely, let us consider some of the masses 
corresponding to internal lines of Fig. la to be large, while the other masses as well as the 
external momentum are small (and of the same order). To distinguish between the small 
and the large masses, we shall denote the latter with capital letters, Mj. By analogy with 
ref. |T7], four different small-threshold configurations exist, namelyQ: 

Case 1: one small 2PT (e.g. masses 7772 and 7775 are small); 

Case 2: two small 2PT's (the masses 7711,1712,1714,171^ are small); 

Case 3: one small 3PT (e.g. masses 1712,1713,1714 are small); 

Case 4: one small 2PT and one small 3PT (e.g. the masses 7772, 7773, 7774, 7775 are small). 
In addition, two special subcases of case 1, when one more mass (not involved in the 
threshold) is small, should be considered separately: 
Case la: case 1 with 7773 being small; 

Case lb: case 1 with one of the masses involved in the second 2PT being small. 

In the present paper we restrict ourselves to cases 1, la, lb and 2 where only small 
two-paiticle thresholds are involved (see a brief discussion of the cases 3 and 4 in the last 
section) . 

Before employing general results on asymptotic expansions in momenta and masses 
| ]T^ , we need to introduce some notation. We shall need Taylor operators Tk and %n 
expanding the denominators in the integrand (in the momentum or the mass, respectively) 

^For the diagram in Fig. lb, pi = p4. If, in addition, mi — 7774, we get the integral with four 
propagators (with one the powers equal to vi + V4). If mi 7^ and both vi and v^i are integer, one can 
use standard partial fractioning to reduce the diagram in Fig. lb to a linear combination of the integrals 
with vi or vi equal to zero. All these integrals can be obtained from the diagram in Fig. la by shrinking 
one of the lines. 

^2PT and 3PT mean two- and three-particle thresholds, respectively. 
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as 

1 ^ (u), (2(M - ey 

i=0 



where 

is the Pochhammer symbol. The "combined" operators T expanding in some momenta 
and masses at the same time can be constructed as products of the above operators (|) 
and (H). It is assumed that the Taylor operators are applied to the integrand before the 
loop integration is performed. After using eqs. (H)-(|D, we should collect together all 
terms carrying the same total power of the small masses and momenta we expand in. 

Let r be the original graph (corresponding in our case to Fig. la), subgraphs of F 
are denoted as 7, and the corresponding "reduced graph" F/7 is obtained from F by 
shrinking the subgraph 7 to a point. Furthermore, is the dimensionally-regularized 
Feynman integral with the denominators corresponding to a graph 7. In particular, Jp 
corresponds to the integral (|^) itself. Then, for our case, the general theorem yields 



k, mj— >0 



where the sum goes over all the subgraphs 7 which (i) contain all the lines with the large 
masses Mj, and (ii) are one-particle irreducible with respect to the light lines (with the 
masses mj). The Taylor operator 7^ g. (see eqs. (|D-(^ ) expands the integrand of 
in small masses mj, the external momentum k and the loop momenta which are 
"external" for a given subgraph 7. The symbol "o" means that the polynomial in g^, 
which appears as a result of applying T to J^, should be inserted in the numerator of the 
integrand of Jr/7- 

The definition of the class of subgraphs 7 involved in the expansion (^ is exactly the 
same as in the zero-threshold expansion |T^, just because zero masses are always small. 
A distinction is that some tadpole-like contributions (vanishing in zero-threshold limit) 
should be taken into account for small non-zero masses (see below). 

Let us consider which subgraphs 7 contribute to the sum (|^) for the different small- 
threshold configurations. To indicate which lines are included in the subgraph, we shall 
use the numbering of lines given in Fig. la. For example, the "full" graph F = {12345} 
includes all five lines, {134} corresponds to a subgraph where the lines 2 and 5 are "omit- 
ted", etc. The subgraphs contributing to the expansion are shown in Fig. 2, where bold 
and narrow lines correspond to the propagators with large and small masses, respectively. 
Dotted lines indicate the lines omitted in the subgraphs 7 (as compared with the "full" 
graph F), i.e. they correspond to the reduced graphs F/7. 

Comparing Fig. 2 with the four first lines of Fig. 3 presented in [|T7| (zero-threshold 
expansion), one can see that Fig. 2 contains some additional subgraphs for the cases la 
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and lb (namely, the second and the fourth subgraphs for each of these cases). These sub- 
graphs did not contribute to the zero-threshold expansion since the corresponding reduced 
graphs r/7 yielded massless tadpole-like integrals and therefore vanished in dimensional 
regularization Here, we get tadpoles with small masses, which do not vanish and 



must be taken into account^. 

Case 1. Two subgraphs contribute to the asymptotic expansion: 

(i)7 = r ^ 



(ii) 7 = {134} 
'd> d"g 



[{k - qY - m2]''2[g2 _ mlY^ 

[{k - pf - MlY^ [{p - g)2 - M|]-3 [p2 _ • ('^^ 

Case la. Four subgraphs contribute to the asymptotic expansion: 

(i)7 = r ^ 



[p2 _ JVf2]^4 

1 



(ii) 7 = {1245} ^ 



d> d"g 



[p2 - mlj'^s 



(iii) 7 = {134} ^ 

1 



d> d"g 



1 

xrfc,,,^3 p _^)2 _ M2]-i[(p - g)2 - mi]-3[p2 _ M|]-4' 



[iv)7 = {14} ^ 

/ / d"» d"(/ ^ 

J J - mlY'i [g2 - mlY^ [{k - qY - mlY^ 



^^fc.P.'? [(p + qf - MlY^[{k -v-qf- ' ^^^^ 

^Appearance of extra subgraphs for the cases la and lb does not contradict the fact that aU the cases 1, 
la and lb correspond to the same small-threshold configuration (one small 2PT). We shall discuss details 
of this connection in section 3. 
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Case lb. Four subgraphs contribute: 

(i)7 = r ^ 



(ii) 7 = {1235} =^ 

1 



j j d> d'^g 



[p^ — ml]''* 



xrfc,p,^,,„, [(p - qf - M|]-3[g2 _ ^2].5p _ p)2 _ M2]-i[(A; - qf - mlY- ' ^^^^ 
(iii) 7 = {134} =^ 

1 



//d>d",. 



1 

[(A; - pY - M2]-i[(p -qf - Mill's [p2 _ ^2].4 ' 

(iv)7 = {13} 



/ / ^""^ 



[p2 - m|]'"4[(/c - g)2 _ mi]''2[g2 _ ^2jl.5 



Case Four subgraphs contribute: 

(i)7 = r 



I l'>' \p-qf-MlY^ 

XT, ^ 



; (16) 



(ii)7 = {134} ^ 

1 



j j d> d"g 



[(A; - gf)2 - mlY^[q^ - mlY^ 



(iii) 7 = {235} contribution can be obtained from the previous one by the permutation 

1 ^ 2,4 ^ 5; 

(iv) 7 = {3} ^ 

//^ ^[(A;-p)2 - mf]-i[(A;-g)2 - ml]-2[p2 _ ^2].4[g2 _ ^2^vs ^'"^ [(p-g)2 - M|]-3 ' 

(18) 



All the formulae (^-(|1^), along with eq. (j^), can be used for any values of the space- 
time dimension n and the powers of propagators i>i. Note that contributions from different 
subgraphs possess ultraviolet and/or infrared poles in e. Cancellation of these poles in 
the sum (when a finite integral is considered) turns out to be a very strong check of the 
calculational procedure^. 

Studying eqs. (|6|)-([T8|) shows that, after partial fractioning is performed wherever nec- 
essary, the following types of contributions can occur in the expressions for the coefficients 
of the small-threshold expansion (in situations with two-particle small thresholds): 

(a) two-loop vacuum diagrams with two (or one) large-mass lines and one (or two) mass- 
less lines; 

(b) products of a one-loop massive diagram (with small masses and external momentum 
k) and a one-loop massive tadpole; 

(c) products of two one-loop massive diagrams with external momentum k. 

The contributions of type (a) are discussed in detail in Appendix A of [O. The 



one-loop two-point functions involving small masses of internal lines (occurring in the 
contributions of types (b) and (c) ) should be calculated exactly, and they are "responsi- 
ble" for describing the two-particle threshold irregularities. Some technical issues related 
to the calculation of these contributions are collected in Appendix A. 



3. Analytic results 

To get analytic results for the terms of the expansion describing the small-threshold 
behaviour, we have used the REDUCE system for analytical calculations The con- 
structed algorithm works for any (integer) powers of the propagators and can be applied 
for both convergent and divergent diagrams. 

As an example, in this section we consider the diagram in Fig. la with unit powers 
of propagators, the so-called "master" diagram. The corresponding scalar integral is 
involved in many interesting physical applications. If all five i>i are equal to one the 
corresponding diagram is finite^ as — > 4, and we shall calculate the corresponding 
results in four dimensions. The algorithm makes it also possible to consider higher terms 
of the expansion in e = |(4 — n), or even results for arbitrary n. 

Let us write the "master" integral a^ 

oo 

J(mi, m2, ms, m^, m^; k) = J(l, 1, 1, 1, 1; mi, m2, ms, 7714, m^; k) = -vr^ ^ Sj , (19) 

j=0 

where 5*^ are the terms corresponding to our expansion. For a given j, the term 5*^ is 
a sum of all contributions of the order {k'^Y" Yl{mfy^ (where the product is taken over 

®See also a discussion in ref. p. 541. 

^To be precise, there is a very special case when this "master" integral becomes infrared-divergent: 
when there is a two-particle zero-threshold and, in addition, — 0. The exact result for this case can 
be easily obtained by a direct calculation. 

"'^'^In eq. dl^), mi may correspond to either small or large masses. In the expressions for the Sj presented 
below, we distinguish the small masses and the large masses by denoting them as rrij and Mj , respectively 
(as in section 2). 
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all lines with small masses) with jo + = j- Since k^J is dimensionless, so should 
be k'^Sj. Therefore, the Sj should involve the corresponding powers of the large masses 
in the denominators. The numerators of Sj may also contain polynomials in the large 
masses, logarithms of the ratios of the masses and some functions corresponding to the 
contributions (a), (b), (c) described at the end of section 2. Note that (in the case when 
there are no zero thresholds) we should not get the terms like ln(— A;^) producing a cut in 
the fc^ plane starting from the origin. 

The contributions of type (a) are expressed through the following dimensionless func- 
tion of two variables: 



n{Ml Ml) = 2U, ( 1 - ^) + 1 In^ ) . (20) 



It is antisymmetric in its arguments and therefore vanishes when Mf = M|. This func- 
tion is connected with the finite part of the two-loop vacuum integral with two general 
masses and one zero mass (for details of this definition, see Appendix A of [|16|)- In the 
contributions considered in this paper the H function depends on the large masses only. 

The contributions of types (b) and (c) are expressed through the function r(mi, k"^) 
which is related to the finite part of the one-loop self-energy with general masses (see 
eq. (|52|)). It is well-known that the r function can be expressed in terms of logarithms 



and square roots (see e.g. in [^). By using the notation 

A(mi, ml, fc^) = Amlml - {k"^ -m\- mlY (21) 
for the "triangle" Kallen function^^, the result for the r function can be presented as 



r(mi, ma; A; ) = 7-^ < V-A In — — — ^= + (m^ - m^) In 



m 



2 



+i7rv^ (k^ - (mi + m^f) | , (22) 



where A = /\{m\,m\,k''). In our contributions, the r function depends on the small 
masses only. In fact, it contains the main information about the small-threshold behaviour 
at two-particle thresholds. The 6 term in the braces yields an imaginary part in the region 
beyond the physical threshold (i.e. for k'^ > (mi -l-ma)^). This is exactly the point where 
the cut in the complex k'^ plane starts. Further properties of one-loop self-energy diagrams 
and the r function are discussed in Appendix A. 

Now, let us present explicit results for some lowest terms of the expansion (0). 

Case 1. For the case of general masses, the Sq contribution to eq. (O) is given by 



So 



4(Mf - Mi) [Mi - Mi){Mf - Mi) 



x<^-4 r(m2,m5;fc^) + |ln \ ^ +2] 



Ml{Ml-Ml) In ^ -Ml{Ml-Ml) In ^ 



^^The function ( pT| ) is totally symmetric with respect to all its arguments. It vanishes at the threshold, 
— (toi + ^2)^: a-iid at the pseudo-threshold, fc^ = (mi — 7712)^. 
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HM'i - M|)(Mi - M|)l„^ + ,n^) - 2M|(M? - M|)ln^ln^ 



+ 2{Mf + Mi){Mi - Mi)n{Mf, Mi) - 2{Mi + Mi){Mf - Mi)n{Mi, M|) I , (23) 
We note that this result is very similar to the corresponding contribution for the zero- 



threshold case, Co (see eq. (21) of ref. jT^). Indeed, in the massless limit we get (cf. 



2 



1-2 




Uln (24) 



Therefore, the massless limit of the brackets in (|23|) containing the r function yields the 
same combination involving ln(— /c^) as in eq. (21) of |T^, while the remaining terms 



coincide. The next term. Si, involves the k , 1112 and mi contributions (we present the 
corresponding result in Appendix B). Note that in this case the coefficient of k"^ gives 
(in the massless limit, i.e. using eq. (0)) the corresponding zero-threshold coefficient 
Ci (see eq. (22) of |T^), while the and m\ contributions have no analogy with the 
zero-threshold case. 

We have also obtained higher terms of the expansion (|T^) with five different masses 
(5*2 and 6*3), but they are more cumbersome, and we do not present them here. 

When some large masses are equal, one can use this fact from the very beginning, 
when calculating contributions of the corresponding subgraphs. Another possibility is to 
consider, with due care (since both denominators and numerators vanish in this limit), 
the corresponding limit of the expressions with different masses. On one hand, the second 
option is a good way to check the consistency of the obtained results. On the other hand, 
the first way (i.e. having fewer different masses from the very beginning) simplifies the 
calculation and therefore makes it possible to obtain expressions for higher terms. 

For example, we get the following expressions for the first two coefficients correspond- 
ing to case 1 with Mi = M4 = M: 

So = - (^2_'m|)^ { (-^-- + ^ 1- (^3 In § + - Ml) 

+Ml (n{M\ Ml) + \ In^ + 2{M^ - M|)| , (25) 



S 1 

' 12M2(M2 - M|)4 

X I fr(m2, m,- e) + \\n ^^"j Um'MI In ^ (k^Mj - {ml + ml){2M^ + M| 
I \ 7712/77.5/ \ M3 



+ {M^-Ml) (k'^{M^-bM'^Ml-2Ml) + ?,{ml + ml)M^{M^ + M 
-6M^Ml I n{M^, Ml) + \^^^J^] {k'^Ml - 2{ml + ml){2M^ + M. 



m^jy^ivi -f M'^ 

+3M^ (-^In^+m^ln^) (2M|(2M^ + M|) In ^-(M^-M|)(M^ + 5M|) 
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+2M^ In ^ (fc2(5M^ - M^) + Vl{m\ + m^^M'^M. 



Ml 



+ (M2 - Mg^) (6(m^ + m^^M^{M^ + 7M|) - k^{M'' + ITM^M^ + 2M 



(26) 



For this case, we have also obtained the 5*2, 5*3 and 5*4 terms. 

If all large masses are equal (i.e. M\ = M3 = M4 = M), the expressions become very 
simple: 



5'n 



2JVP 



r(m2,m5; /c ) + ^ In 



2 2 



+ 3 



(27) 



^1 



1 



144M4 



6(A;^ + ml + ml) r(m2, ms; /c^) + | In 



2 2 



+6 ( m2 In— J + ml In — ^ 1 + 9r + 26(m^ + "^5) | • 



(2^ 



For this case, we have obtained the terms up to (and including) 5*6. 

Case la. For general masses, the lowest case la contributions to (|T^) are 



Ml - Ml Ml 



. Ml ( , , 2x 11 Ml Ml 

In ^ I r(m2, mg; fc^) + ^ In !, t + 2 



2 2 

7722^5 



(29) 



^1 



4M2M|(M2 -M|)3 



2 r(m2,m5; A;^) + Mn 



2 2~ 

m2m5 



+ 2 



M^Ml ({Ml + M2)A;2 + (M^ - Ml){ml - m^)) In - 2{Ml - Ml)k 



Mi 

HMl - Mlf Uml [mi In ^ - Ml In ^1 - + Mlml 



ml 



M I m 
^MlMUMl - Ml) 1„ ^ ™^ in ^ - 1„ 



- 2m2 + 2ml 



■{Ml -Ml 



Mlml In 
Ml 



mT, 



Mlml In ■ 



ml 



MlMl ^'^^ '"^MlMl 



+ {ml + ml){Ml + Ml)\n 



Ml 
Ml 



+2mi In ^ ^M/ In ^ + M,^ In ^ j + 4..^ ^M/ In 
+ (M2 - Ml) (A;2 - m: 



- In 



+ 4m^(|7r2 - 1) ) - 2M^m^ + 2M^m^ 



.(30) 



For the general case la, we have also obtained the terms up to (and including) 5*5. 
If Ml = M4 = M, the corresponding results are 



5n 



1 



r(m2,m5; fc^) + 5 In 



M^ 



2 2 
m^m 



+ 2 



(31) 
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^1 



12M4 



r{m2, ms; k'^) + 5 In —232 J ( +3(m^+m^+4m|) + 12m| In 



2 2 



i 

M2 



o / 2i "^2 2i 

-3 I TTio In ^ + 771^ In ■ ^ 



M2 



+ 6(771^ + 777^ + 4777^) - 47r^777^ } .(32) 



Case lb. The results for the lowest terms, So and Si, are given by 



2Mf (Mf - M|) 



2M2 In ( t(7772, 7775; A;') + i In 



77727775 



+ 2 



2 M|_ 



-|7r^(M,^-M3^) ,(33) 



Si = 



4MtMi{Mf - Mif 



-2Ml (t(7772,7775;A;^) + |ln 



+ 2 



X 



/ ^#2 

M.^Ml ( {{Ml - Ml){e - ml) + {Ml + )777g In ^ - 2{Ml - Ml)m 
+ {Ml - Mlf ( 2777^ I M^n^ - Mfm^) " " " 



Mi 



Mie + 



+M| [n{Ml,Ml 



,,^2Mr 

2 



(Mf - M^Y{Ml + M3^)(r + 2777^ 



+ 2MlMl{Ml - Ml)k'^ + AMlMl{ml + 777^) 



+{Ml-M'i)k'' i2MlM'i{Ml + M'i) In ^ + Mi^(Mi'-3M|) + \'K''Ml{Ml-Mlf 



+A{Ml - Mlf ml I In ^ - Ml + W{2Ml + M, 



777^ 



Mi 



+2Mt -{Ml - Ml)ml In ^ M| In ^ - Ml + M| 



Ml 



+ Mo 777^ In 



Ml 



mi 
Ml 



Mi 

' M^ 
{Ml^Ml)\^-l^-2Ml 



mr 



2M' 



+2MfMi In [{2Ml + Mi)mi + Mi{k' 

+M^{Ml - Ml) ({hMl - 3Ml)ml + {7 Ml - MDmj) | . 

For the general case lb, we have also obtained the S2, Ss and S4 terms. 
If Ml = M3 = M, the results for Sq and are 

M^ ^ 



'S'o = --^ (r{m2, 7775; A;^) + | In 



2^2 + + 4 



7775777 



(34) 



(35) 



Si^ 



12M4 



r(7772,7775; A;^) + |ln 



M^ 



2 2 



+ 1 



3k'^+m^+3ml+12ml In 



777^ 



M2 



+ 1 



-3777^ In ^ - m^ln^ - A;^(7r^-10) - 6777^(7r^-8) + ^(777^+777^) | .(36) 
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Before proceeding to the case 2, let us discuss a connection between the results for 
the cases 1 and la,b. All of them correspond to the same small-threshold configuration 
(one small two-particle threshold), and one could expect that there should be a transition 
between explicit expressions for the terms of the expansion ([T9| ) for the case 1 and those 



for the cases la and lb. In the zero-threshold case [|17|, what we needed was just to put 
one more mass (in the coefficients Cj for the case 1) to be zero, and we arrived at the 
results for the corresponding case la and case lb coefficients. When we deal with small 
(but non-zero) thresholds, the situation is more tricky. We need to consider one more 
mass to be small, but now the contributions proportional to this small mass will go to the 
higher Sj. In addition, we need to expand the denominators involving this small mass as 
well as the corresponding Ti. function(s) involving a small and a large mass. Expansion 
of the denominators is trivial, while the Ti function can be expanded using the following 
formula: 

So, if one wants to show the correspondence between the case 1 and the cases la,b 
, one should not consider the Sj themselves but to collect the terms of required order 
in the whole sum (|19D instead. In such a way, we have checked that all available terms 
for the case 1 produce the corresponding terms for the cases la and lb. On one hand, 
the described additional "subexpansion" in a small mass corresponds, in some sense, to 
additional subgraphs appearing in the cases la and lb as compared with the case 1, see 
Fig. 2. On the other hand, if one is interested in cases la and lb, it is better to consider 
these cases from the very beginning rather than get results from the case 1 (since more 
terms are available for the cases la and lb). 

Case 2. The results for the 5*0 and 5*1 terms of the expansion ([T9|) are 

•^0 = -i;72 I r ("^1 , "^4; A;') + i In + 1 ) ( r(m2, m^-k'') + \\ii + 1 ) + ivr' - 1 

M3 [ y ^ rriimi J \ ^ 7713^5 / 



o ^ ( n 2 2 2 2 2\ 



r(mi,m4; /c^) + iln^^^-| ) ( r(m2, ms; fc^) + | In ^ 




9 2 2 



r(mi, 7714; A;^) + | In — - || ( In ^ + In ^ - |(m^ + m^) + A;^ 



+ In ^ - |(m2 + ml) + k'^ 



0)) (r(m2, mg; /c^) -r(m2, mg; 
where the value of the r function at = is 




mi,m2;0 = -1 - 5 ^ 40 

^ mf — 7712 ^2 



13 



For the general case 2, we have obtained the terms of the expansion (|T9D up to (and 
including) Sq. 

An interesting feature of eq. (|39D is the appearance of k"^ in the denominator. For 
higher terms, higher powers of k"^ in the denominator occur^ e.g. (fc^)^ for S2, etc. For 
general masses, the powers of fc^ in the denominator can be cancelled by considering the 
Taylor expansion of the r function in fc^. The simplest example can be seen in eq. (^9]). 

To conclude this section, we would like to note that we have successfully compared 
the massless limit of the small-threshold expansion with all coefficients from the four first 
columns (corresponding to the cases 1,1a, lb and 2, respectively) of Table 1 presented on 
p. 545 of ref. 0. 



4. Numerical results 

Let us show how the small-threshold expansion can be applied to obtain approximate 
numerical results for self-energy diagrams. We shall use as examples a number of mass 
configurations corresponding to diagrams occurring in the Standard Model. We take the 
masses of the W boson, and of the charmed, bottom and top quarks to be[^ 

Mw = 80 GeV, = 1.5 GeV, mi, = 5 GeV, = 174 GeV. (41) 

In this section, we shall consider the following approximations to the "master" integral 
(P^D (corresponding to Fig. la): 

j(N) ^ j2 S, . (42) 

j=0 



As in ref. IjT^, the first example corresponds to a diagram containing a top-bottom 
loop to which two W bosons are attached. This diagram contributes to the self-energy of 
the photon and the Z boson. The corresponding scalar integral is 

J{Mw,mb,mt,Mw,mb;k) . (43) 

If we consider rrify as a small mass, and rrit and Mw as large massesQ, this diagram has 
one small two-particle threshold at k'^ = 4ml therefore it belongs to case 1. In Fig. 3 
the approximations J^^^ defined by eq. (^2|) are shown as curves, and, for comparison, 
values of J obtained by numerical integration are shown as crosses. The position of 
the lowest "large" threshold, at k^ = 4:M^^ in this example, is indicated by a vertical line. 

"'^^The corresponding contributions vanish when either = m| or = mg. 

-'^^In this section, we adopt "standard" notation for the masses of physical particles: the quark masses 
are denoted with small m whereas the capital M is used for the vector boson masses. To avoid any 
confusion with the notation of sections 2 and 3, we shall explicitly state, for each concrete example, 
which masses are considered to be small and large. 

^'^Note that in ref. [l^ we have considered the situation when mi, = 0. The situation for non-zero mi, 
was discussed as a subject to be investigated in the future. This paper solves the problem posed in [p^. 
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One can see that the expansion converges all the way up to the first large threshold at 
A;2 = 25600 GeV^. 

In fact, the behaviour for Ami ^ < 4M^ is described as well as in ||T^. The main 
difference is in the region of small k"^ (in particular, around the 4m^ threshold) which 
could not be described by the zero-threshold expansion [0. In the scale of Fig. 3, the 



details of the small-threshold behaviour cannot be seen so that we use a different pair of 
plots in Fig. 4 to illustrate the behaviour in the region of small momenta, including the 
threshold at fc^ = 4mg = 100 GeVl One can see that even the zero order approximation 

J(0) 

is very good and that it is sufficient to take J^^^ to get result with a high precision. 
We do not present other approximations (up to = 4) because they give a precision 
better than the numerical program does (and one cannot distinguish them in the plot). 
In our approximations, the irregularities at the small threshold are reproduced by the r 



function (22). In particular, the imaginary part is zero below the threshold and behaves 
like a square root in the region above the threshold. 

In fact, one can get a better description of the behaviour of the same diagram (|43|) in 
the region around the fc^ = 4M^ threshold by considering as a small mass. In this 
case, the only remaining large mass is rrit. Thus, it belongs to case 2, in our terminology. 
According to our results, the situation in the region of small momenta is more or less the 
same: the first two approximations happen to be good enough (but not so good as when 
the diagram is considered as case 1). Our approximations in the region up to the large 
three-particle threshold (at k'^ = (m^ -|- M^y + m^)^) are shown in Fig. 5. One can see that 
our approximations describe the behaviour around the 4M^ threshold and, moreover, 
they work even beyond this second threshold^. Thus, just by treating in another way the 
available massive parameters it is possible to characterize analytically the first diagram 
in a larger region of momenta. 

As second example we consider the integral 

J{mb,mc,Mw,mb.mc;k) (44) 

which has two small two-particle thresholds, at k"^ = 4m^ and k"^ = 4m^, and therefore 
belongs to case 2. 

The behaviour up to the first large threshold is rather similar to the previous examples 
so we do not present the corresponding plot. The behaviour at small momenta including 
both small thresholds can be seen in Fig. 6. A rather tricky behaviour in the region 
around these two small thresholds is very well reproduced by taking the lowest analytical 
approximations, J*^"-* (for the imaginary part) or J^^^ (for the real part). 

All the examples illustrate that the small-threshold expansion provides approximations 
which perfectly describe the small-threshold behaviour. Moreover, they are at least as 
accurate as numerical integration in a large part of the region of convergence and can be 
evaluated much faster. 



^^The convergence between the 4M^ and the (mt + AIw + rnt)'^ thresholds is not so good since the 
Mw I rrit ratio is not very small. 
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5. Conclusions 



In this paper we have studied the behaviour of two-loop self-energy diagrams in the 
situation when the external momentum and some of the masses are small with respect 
to the large masses. All configurations with small two-particle thresholds have been 
considered. By use of explicit formulae for the terms of asymptotic expansions in the 
large mass limit, we presented an analytic approach to calculating these diagrams by 
keeping the first few terms of the expansion. 

By taking some complicated cases (corresponding to diagrams with different masses 
occurring in the Standard Model) as examples, we compared our results with those of a 
numerical integration program based on the algorithm of ref. [0 (see also in |^). We have 
shown that our analytical approximations work very well in the region of small momenta 
that includes all the small thresholds. Moreover, the small-threshold expansion converges 
up to a region close to the first large threshold. Unless fc^ is close to the large threshold, 
only a few terms are needed to obtain accurate results. This comparison can also be 
considered check of the numerical program. 

Thus we have solved the problem of the threshold behaviour of two-loop two-point 
functions for the small two-particle thresholds. It is interesting that it was possible to do 
this just by using the expansion in the large masses. The main idea was to avoid putting 
any conditions on relative values of the external momentum squared and small masses. 
To our knowledge, no efficient algorithms for describing the non-zero threshold behaviour 
of two-loop diagrams with arbitrary masses were available so far^ A nice feature of the 
presented approach is that all two-particle threshold irregularities of two-loop diagrams 
are incorporated in the function corresponding to one-loop diagrams, the r function. 

A natural extension of this work could be connected with cases 3 and 4 (according to 
the classification given in Section 2), when three-particle small thresholds arise. To carry 
out a similar program for those cases we can apply the same technique of asymptotic 
expansion in large mass(es). In particular, in eq. (^) one should sum over the subgraphs 
obeying the same conditions as in the cases considered. For example, for case 3 (with 
small masses 7712,7713 and 7714) the following five subgraphs 7 contribute to the sum (|^): 
r = {12345}, {1235}, {1345}, {1245} and {15}. For case 4 (when all the masses except 
Ml are small), four subgraphs contribute: F, {1245}, {134} and {1}. However, in these 
cases the small-threshold behaviour at three-particle threshold is defined by the functions 
corresponding to the sunset diagram (with three propagators) and the diagram with 
four propagators (e.g. the integral (^ with ui = 0). Unfortunately, sufficient analytic 
information about threshold behaviour of these diagrams is at the moment not available. 

The algorithm presented in this paper can be also extended to the three-point two-loop 
diagrams with different masses. 

^^A possible exception may be related to an interesting approach considered in ref. ||2^ where a two- 
loop diagram with one non-zero mass parameter was studied as an example. It is not clear, however, how 
efficiently the method of will work for two-loop diagrams with different non-zero masses, since the 
coefficients of the expansion may require results for the threshold values of two-loop integrals (which, for 
the general case, are not known analytically). 
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Appendix A. One-loop two-point integrals with masses 

In this appendix, we collect some relevant results for the integrals 

Let us introduce dimensionless quantities 



mi 



ms 



X{x,y) = X - yY - Axy. 

It is easy to see that the Kallen function A, eq. ([21|), is related to A via 

A(mi, ml, k^) = Amlml - {k'^ - ml- mlf = -{k'^f X^{x, y). 



(45) 

(46) 
(47) 

(48) 



Therefore, A also vanishes at the threshold, k"^ = (mi+m2)^, and at the pseudo-threshold 
fc^ = (mi — m2)^. 

For arbitrary z/i,z/2 and the space-time dimension n, the result for the integral 
can be written as (see in [Q) 

J(z/i,z/2;mi,m2) = vr"/¥-"(A;2)"/2~-i-2 

r(t-^i)r(t-z.2)r(z.i + z.2-t 



r(z/i) r(z/2) r(n - i^i- V2) 

XF4 (z/i + z/2 - f , i/i + 1/2 - + 1; z/i - f + 1, z/2 - f + 1 

r [v, - 1) 



+ (-x 



r(^i) 
r(^2) 



(^^2, - f + 1; f - + 1, - f + 1 

z/i - I + 1; z/i - I + 1, I - zy2 + 1 



x,y 
x,y 

x,y 



(49) 



where F4 is the Appell hypergeometric function of two variables, 
and {a)j = r(a + j)/r(a) is the Pochhammer symbol. 



(50) 
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For z/i = z/2 = 1, but keeping the space-time dimension n = 4 — 28 arbitrary, the result 
simphfies and can be written in terms of the Gauss hypergeometric function 2-^1, 



2-e I 7,2 



J(l, 1; mi, = i VT {~k 



■Qi^£)r(£) 

r(2 - 2e) 



-2e 



-r(-i + i(i + X - y - \)i-x)~' 

-r(-l +e)l{l-x + y- X){-y)-' 2^1 
Expanding in e = | (4 — n) we get 



1,6 
2-6 

1,6 
2-6 



{l + x-y-Xy^ 
Ax 

(1 -x + y- Xf^ 
4|/ 



■ (51) 



J(l, 1; mi, m2) = ivr^"" r(l + |^ + 2 - i (Inm? + lnm|) + r(mi, m2; k^) + . 



where the r function can be defined via an integral representation 



r(mi, m2; fc^) = r(x, y) = — da \n 



It is understood that ^ + iO. 



' am\ + (1 — a)m\ — a{\ — a)k'^ 
a(l — a;)mim2 



(52) 



(53) 



For fc^ < (mi — m2)^ (including the Euclidean region, k"^ < 0), the result for the r 



function can be written as 



^(a;, y) = |A In ] — - — ^ + |(x - y) In -. 

I — X — y + X ^ X 



(54) 



This is also valid all the way up to the threshold, fc^ = (mi + m2)^. However, since X{x, y) 
is imaginary for (mi — m2)^ < k"^ < (mi + m2)^, < 0, in this region one can also use 
another representation. 



arccos 



' X + y — V 



Ux-y) In- 



X 



(55) 



Beyond the threshold, at A;^ > (mi + m2)^, is again positive, but the r function obtains 
an imaginary part, 

r(x,?/) = ivrA + iA In^^ — - — - — ^ + l(x - ?/) In -. (56) 
^ 1 — X — y + X ^ X 



Note that at A = the expression for the r function is 



T{x,y) 



\{x - y) In-, 



A=0 



X 



(57) 



whereas at fc^ = we get the result (^0|). 

When one or two masses m^ vanish, the r function develops a logarithmic singularity. 
For example, taking the limit m2 — >■ yields the well-known result 



r(mi,m2; fc^ 



m2— »0 



1 , ml rn^ — k"^ , — k"^ 
lln — In — . 

^ m? k"^ nii 



(58) 
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The term | Inml cancels the corresponding term in (|52D. The rest is real for k"^ < ml, i.e. 
below the threshold (in this case the threshold and the pseudo-threshold coincide). For 
k"^ > ml (beyond the threshold), we should remember that k"^ k"^ + iO and substitute 
ln(mf — fc^) <-> ln(/c^ — mj) — ivr. 
If both masses vanish, we get 



1 A. f "^'A , 1. ^ ^2 



2 

mi ,m2—'0 



The terms |(lnm^ + Inmg) cancel the corresponding terms in (|5^), while the remaining 
ln(— /c^) should be understood as (InA;^ — ivr) for the time-like values of the momentum. 
In this limit, the results of the present paper correspond to the results of zero-threshold 
expansion |jl7|. 

The integrals J(z/i, z/2; mi, 777,2) with higher integer values of z/i or z/2 can be obtained 
by using the integration- by-parts technique p5|. The way is quite similar to one used in 
[p6l . First, we get a system of two equations, 

- 2i'im\j{i'i + 1, z/2) + z/2(A;^ -m\- m^J^Ui, V2 + 1) 

= -{n - 2z/i - U2)J{ui, z/2) + z/2J(z/i - 1, z/2 + 1), 

iyi{k'^ ~m\- m\)J{vi + 1, z/2) - 2v2ni\j{vi, z/2 + 1) 

= -{n - z/i - 2Z/2) J(z/i, z/2) + z/i J(z/i + 1, z/2 - 1). (60) 

The determinant of the matrix composed of the coefficients on the l.h.s. of this system is 
proportional to A{m1,ml, k"^), eq. (^Sj). Solving the system, we get 

J(z/i + 1, z/2) = ~ "^i)(^ - z/i - 2Z/2) + ml{n - 3z/i)) J(z/i, z/2) 

-2z/2m^ J(z/i - 1, z/2 + 1) - z/i(/c2 - - m^) J(z/i + 1, z/2 - 1)} , (61) 



and an analogous result for J(z/i,z/2 + 1). By using these results, all the integrals with 
higher integer z/'s can be expressed in terms of the integral J(l, 1; mi, m2) and the massive 
tadpoles. 

We also need some formulae for the integrals with numerators. For the integrals ( ^3D 
with one negative power of the denominator, we get 

J{iy,-N;mi,m2) = J ^^2 1 1^2^ [*^^ ~ ~ "^2] ^ 

Another formula we would like to present is useful for dealing with the two-loop 
integrals with the denominators corresponding to a product of two one-loop integrals (|^), 
but when the numerator is a power of the scalar product of two integration momenta. 
We need such integrals to calculate the terms of the expansion (|19D for case 2. We shall 
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present the result valid for a more general case, when the integrand is a product of two 
arbitrary scalar functions. Let 

[something] = J jdTpdTq {something} Fi{k'^ , p"^ , (pk)) F2{k'^ , , (qk)) . (63) 



Then 



K 



(pq) 



N 



(!) 



X K 



(64) 



where the sum goes over all {j, js} such that 2j + = N (so, in fact this is a one-fold 
finite sum) while^'^ 



1=0 



l\{j-2l)\ 



(65) 



are Gegenbauer polynomials. The result ( p^ ) can be derived using the formulae presented 
in ref. ||2^|3- Using the definition (|65D , it is easy to see that we do not get neither negative 
powers nor square roots of and on the r.h.s. of eq. 



Appendix B. The 5*1 term for the case 1 with general masses 

Here we present the result for the Si contribution to the expansion ([T^) for the case 1 
with different masses]^: 

1 f / M^M^ \ 



X 



MliMl - Mlf In ^ {MliMl - Mi){Ml - Ml){k^ + - ml) 

+2Ml{Ml -Mlfk'^ + 2Ml{Ml - Mlfml 

+Ml{Ml - MD^ln ^ {Ml{Ml - MI){MI - Ml){e -ml + ml) 

-2Ml{Ml - Ml)^k^ - 2Ml{Ml - M^)^ml 
-{Ml - Ml){Ml - Ml){Ml - Ml) 

X ({Ml - Ml){Ml - Ml) (mI{MI - Ml) + Ml{Ml - Ml 



^^The symbol [j/2] denotes the integer part of j /2. 

^^li is interesting to note that eq. ( p^ has a structure similar to eq. (43) of with Ni — N2 = N 
and ji = j2 = j- 

-^^The lowest term, Sq, is given by eq. (p3[). 
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+ {Ml - Ml){Ml - Ml) (mI{MI - Ml) - Ml{Ml - Ml)) mj 
-{Ml - Ml){Ml - Ml) {mI{MI - Ml) + Ml{Ml - Ml)) mfj 
^ M^\ r 

HiMl, Ml) - I In' j^j {Ml - Ml)' [4MlMl{Ml - Mlf {ml + ml) 

^Ml - Ml) ({Ml + Ml){M^ - Ml) + 2MlMl{Ml - Ml)) k'^ 



m: 



+ [n{Ml, Ml) - \ In^ ^ j {Ml - Mlf [4MlMl{Ml - Mlf {ml + ml) 

+{Ml - Ml) [{Ml + Ml){Ml - Ml,) - 2MlMl{Ml - Ml)) k' 



+2{Ml - Ml){Ml Ml)ml In ^ 



M2 



{Ml-Mlf M^ In {Ml -Mlf (mI{MI-MI)-2MI{MI-MI)) Ml\n 

- {Ml - Ml){Ml - Ml){Ml - Ml) (mI{MI - Ml) + 2MI{MI - Ml 

+2{Ml-Ml){Ml-Ml)ml]n^ 
M' 

{Ml-Mlf Ml In j^- {Ml -Ml)'' (mI {Ml -Ml) + 2 Ml {Ml -Ml)) Mlln 
{Ml - Ml){Ml - Ml){Ml - Ml) [mI{MI - Ml) - 2MI{MI - Ml 



Ml 
Ml 



X 



Ml 



Ml 



+ 



-Ml {Ml - Mlf In' ^ \Ml{Ml - MI){MI - Ml){k' + - ml) 



+2Ml{Ml ~ Mlfk' + 2Ml{Ml - MD'mt 



-Ml{Ml - Mlf In' ^ [Ml{Ml - MI){MI - Ml){k' -ml + m',) 

-2Ml{Ml - Ml)'k' ~ 2Ml{Ml - Mlfn^ 
-Ml{Ml - Mlf In ^ In ^ [mI{MI - MI){MI - Ml){k' -ml- m',) 

-2Ml{Ml - Mlfml - 2Ml{Ml - MD'ml 

M^ r 

+{Ml-Ml){Ml-Ml) In -I [Ml{Ml - Ml){k' -ml + m',) 



X 



Ml - Ml) (2Ml{Ml - Ml) + Ml - M^) - 2{Ml + MI){MI - Mlf 



II? 



-{Ml - Ml)' (8MlMl{Ml - Ml) + 3Ml{Ml - M^) - Ml{Ml - M. 
+{Ml - Ml) (2Ml {Ml - Ml)'{Ml + 2Ml) - {Ml -Mlf {Ml - MI){MI+MI)) ml 

M' r 

+ {Ml -Ml) {Ml -Ml) In -j^ [mI{MI - Ml){k' + ml - ml,) 

X ({Ml - Ml) (2Ml{Ml - Ml) - Ml + M|) - 2{Ml + MI){MI - Ml)') 
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+ {Ml - Mlf (8MlM^{M^ - Ml) - 3Ml{M^ - M^) + Ml^Mj - M. 
+ {Ml -Ml) {2 Ml ( Ml -Mlf{Ml + 2Ml) + {Ml-Mlf{ Ml -Ml){ Ml + M. 
- (Ml -Ml) {Ml -Ml) {Ml -Ml)"^ [{Ml - MI){MI + Mf)(M| + Ml){k'^ -ml-ml) 

-AMI {Ml -Ml) {Mlml+Mlml + Ml {k'^- 2ml -2ml 

-2{Ml - Ml) {Ml - Ml) {Mlml - Mlml) 



.(66) 
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Figure 1: Two-loop self-energy diagrams. 



Case 1 
Case la 
Case lb 
Case 2 




-Q 






-Q - < 





Figure 2: The subgraphs 7 to be included in the sum (5) 
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Figure 3: The real and imaginary parts of J(M^, mfc, mt, M^, mf,; fc). 
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Figure 4: The real and imaginary parts of J(M^, mfc, mt, M^, mf,; fc) at small momenta. 
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Figure 5: The real and imaginary part of J(M^, mfc, mt, M^, mf,; fc) up to the three- 
particle threshold. 
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Figure 6: The real and imaginary part of J[mb,mc, Mw,'mb,'fnc',k) at small momenta. 
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